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Abstract 

In this work we explicitly calculate the thermal conductivity for a general bidi¬ 
mensional dilute gas of neutral molecules by solving Boltzmann’s equation. Chapman- 
Enskog’s method is used in order to analytically obtain this transport coefficient to 
first approximation in terms of a collision integral for an unspecified molecular inter¬ 
action model. In the particular case of hard disks interactions, the result is shown 
to be consistent with previous work by J. V. Sengers [1], yielding the expected 
dependence with the temperature. This dependence is widely used for dense gases 
as the low density limit in the Enskog expansion. The case of Maxwellian molecules 
is also explored where a linear dependence with the temperature is obtained. 
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1 Introduction 

Most theoretical work on the field of bidimensional gaseous systems is limited to the 
case of dense gases. Moreover, the question of whether a density expansion is appro¬ 
priate and the origin of the divergence of the corresponding transport coefficients 
for low density limits have been topics of intense research [31H1 m I [31 ig. 
To the authors’ knowledge, the only work in which transport coefficients for the 
dilute bidimensional gas are calculated can be found in Ref. [T]. In that work, only 
the particular case of a hard disks gas is presented as a side calculation in order to 
address the triple collision contribution in such a system. 

In the present work an explicit calculation of the thermal conductivity for a 
dilute bidimensional gas is developed in detail without specifying a particular inter¬ 
action potential. The heat flux and corresponding conductivity are then obtained 
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for two particular molecular models, hard disks and Maxwellian particles, the for¬ 
mer yielding results consistent with Ref. [T]. It is worthwile to point out that these 
resultas cannot be directly deduced as a particular case of the 3D one. In particular, 
the dependence of the scattering cross-section with the scattering angle cannot be 
obtained as a reduction of the corresponding in three dimensions. These expressions 
can be introduced in bidimensional numerical simulations and are feasible of gener¬ 
alizations, for example to the relativistic scenario. Additionally, the formalism here 
presented in detail, can be useful in the calculation of other transport coefficients. 

The rest of this work is organized as follows. In Section 2 we briefly review the 
basic formalism for calculating the heat flux from the linearized Boltzmann equation 
within the Chapman-Enskog approximation, here in the case of a two dimensional 
gas. The general expressions for the heat flux and thermal conductivity are obtained 
in Sect. 3 in terms of a collision integral which is calculated to some detail in the 
Appendix. Section 4 addresses the particular cases of hard disks and Maxwellian 
molecules. The final section includes a brief discussion and some final remarks. 


2 Boltzmann equation and Chapman-Enskog approximation 


The evolution equation for the single particle distribution function, corresponding 
to a dilute bidimensional gas of neutral non-interacting molecules is given by the 
Boltzmann equation [lOj 

^ + v-Vf = J [ff[ - ffi] gadxd^vi. (1) 


Here /, the distribution function, corresponds to the particle density in the four 
dimensional phase space such that f {r,v,t) (PrcPv yields the number of particles 
at time t with a position vector between r and f+ dr and velocity between v and 
V + dv. The right hand side of Eq. © accounts for changes in the single particle 
distribution function due to collisions. Non-primed and primed quantities denote 
values before and after the collision respectively. The molecule’s relative velocity is 
denoted by g, with g being its magnitude. Also x a are the 2D counterparts to 
the solid angle and the scattering cross section in the 3D case respectively. However 
it is worthwhile to mention at this point that a (which has dimensions of length) 
depends on g and x even for a hard disks model. 

The local equilibrium solution of the Boltzmann equation for a bidimensional 
system can be shown to be given by 


/( 0 ) = 


nm 

2TTkT^ 


( 2 ) 


by means of which one can establish the corresponding expressions for the state 
variables, in this case the particle density 


the hydrodynamic velocity 


n = 


nu = 


fd?v, 

(3) 

fvd?v, 

(4) 
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and the internal energy density 



The relationship between the molecular (v), hydrodynamic (u) and peculiar (c) 
velocities appearing in the previous expressions is given by T = u + c being v the 
velocity of a molecule as measured by an arbitrary observer while c is the one seen 
from a system comoving with the fluid, that is an observer moving with velocity u. 

At this point it is convenient to introduce the Chapman-Enskog expansion. Such 
a procedure allows one to express the complete solution to Eq. m as an infinite series 
of terms ordered in the Knudsen parameter. This parameter quantifies the deviation 
from the local equilibrium distribution function due to the gradients present in the 
system [10]. To first order in the gradients, corresponding to the Navier-Stokes 
regime, one thus has 

/ = /(°) (! + (/.), ( 6 ) 

In this context, as is well known, the balance equations for the variables given in Eqs. 

can be readily established by multiplying Eq. o by a collisional invariant 
and integrating over velocity space. Such a set requires closure, or constitutive, 
relations expressing the dissipative fluxes in terms of gradients of the state variables. 
In particular, for the heat flux one obtains 

9 = ycd^c, (7) 

which will be used in this work to obtain an analytic expression for the thermal 
conductivity k from the Fourier relation q = —AVT. 

The solution for cj) is obtained from the linearized Boltzmann equation to first 
order in Knudsen’s parameter which is given by nni 

where the linearized collision kernel reads 

= j Wi + (j)' - (1)1 - (t>] gadxdvl- (9) 

Following the usual steps in the Chapman-Enskog method uni, the vector part of 
Eq. ® can be shown to be given in this case by 



where we have neglected the terms involving the gradient of the hydrodynamic 
velocity since it will not be coupled with the heat flux. 

3 Solution to the integral equation 

The Chapman-Enskog hypothesis together with the linearized Boltzmann equation 
yield the integral equation (1101) . Although it is much simpler to solve than the 
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integro-differential equation m, its complete solution can be involved. However, as 
will be shown, the two-dimensional calculation is straightforward and the collisional 
bracket corresponding to the heat flux can be calculated directly for a hard disks 
collisional model. 

The solution to Eq. OT, following the same procedure as in the 3D case, can 
be written as 

(j) = A{c)c-'^, ( 11 ) 

once the subsidiary conditions are enforced m- The function .4(c) needs to be 
determined such that it satisfies the inhomogeneous equation (nni, being still subject 
to the following condition, which also arises from the subsidiary conditions, 

J f^°^A{c)ccd^c=^^-I j /(°U(c)c^(i2c = 0, (12) 

where I is the identity matrix. Introducing Eq. dm in Eq. ra and defining an 
auxiliary variable x = mc^/2kT , the integral equation for A (x) can be written as 


/(O) _ 2) c. ^ = I (^A {x) c • 

(13) 

where A (x) satisfies 


J A{x)p (x) dx = 0, 

(14) 

with 


p (x) = xe~^, 

(15) 


being the weight function. Following the standard procedure, A {x) is written as an 
infinite series in orthogonal polynomials. In this case, as in the three dimensional 

(k) 

one, the appropriate basis are the Sonine polynomials £„ . Taking into account 
Eqs. (fT!?l) - (fT5l) . the appropriate polynomials for the 2D case correspond to A: = 1 
since 

J (x) (j.) (16) 


with 


^0 ( 2 ^) = 1 £[ (x) = 2 — X. 

(17) 

Thus, we have 


A(x) = J2 ' ( 2 ;) > 

n—0 

(18) 

and the condition 

in Eq. (1141) is written as 



CXD „ 

^a„ p {x) (a;) dx = 0, 

n=0 

(19) 


which implies oq = 0. Thus, the solution to the integral equation m is given by 
Eq. (ITT]) with A (x) = (x) and a„ obtained from 


OO 

- {x)c=Y^ anX {x) c) . (20) 

n=l 
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With the proposed solution, the heat flux can be written as 


q = — [kTf — ^ a„ / C^^'> (x) xp (x) dx. (21) 

m 1 I 

n=l 

Writing x in terms of the Sonine polynomials as a; = 2 £q^^ (x) — (x), and using 

the orthogonality condition (HSl) one obtains that 


q = -2-(kTf^a^. ( 22 ) 

m J 

Thus, only oi is required for the purposes of this work, the scalar product of Eq. 
(1^ with d'm (x) c is calculated and the resulting equation is integrated over velocity 
space. The left hand side yields 

f {x)Cimx)p{x)dx = (23) 

m J m 

and thus 

= - ^ a„ f (x) c • I (x) d^c. (24) 

Equation (1241) can now be solved for oi to different orders of approximation using 
the standard variational method, whose details can be found in Ref. m- To first 
approximation one obtains 



The integral above is usually written in terms of a collisional bracket and requires, in 
most cases, tedious calculations in order to be evaluated. However, using momentum 
conservation, one can show that 



(26) 


and by interchanging particles after and before the collision as well as dummy indices, 


j H-I (G'^d‘^c= j G-I 
Using Eqs. (1^ and one obtains 


(-) 


d^c. 


(27) 


where 



and thus oi can be expressed as 


oi = —4- 


nk'^T^ 


r-l 


(28) 

(29) 

(30) 
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The detailed calculation of I is shown in the Appendix where the expression 


I 


n^/3 



g^sin^ Xa- (x,g) dx(fg, 


( 31 ) 


is obtained with /3 = m/2kT. It is worthwhile to point out that in the bidimensional 
case, this integral is in principle simpler to calculate than in the 3D case. The general 
expression for the thermal conductivity of a dilute bidimensional gas is thus given 
by 


A = - 


1 8n2fc4r3 


(32) 


This expression for the transport coefficient is the main result of this work. In the 
following section, two particular models for the intermolecular interaction will be 
introduced in order to obtain explicit temperature dependence of A. 


4 The hard disks and Maxwellian molecules models 


In order to perform the integration in Eq. m with respect to the scattering angle 
X, a particular model for the 2D scattering cross section needs to be introduced. For 
example, in the case of hard disks a (x, g) depends on the scattering angle as follows 



(33) 


where b is the impact parameter, this leads to the following expressions for I and 
oi (see Appendix): 


I = 





(34) 


and 


1 / m 

^ nd \ irkT J ’ 

respectively. Using the results above, the heat flux can be written as 


(35) 


q = 


2 

d 




(36) 


from which one obtains the thermal conductivity for the hard disk dilute gas as 




(37) 


consistently with the results quoted in Ref. [T] . 

The Maxwellian disks model is a particular case of the point centers of repulsion 
potential and approximates the molecular interaction for the case of high tempera¬ 
ture gases in which the repulsive forces dominate over the attractive ones. In such 
a case, the potential can be expressed as m 

(j) {r) = (38) 

u — 1 
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where z/ is a measure of the hardness of the molecules. Maxwellian molecules are 
characterized by ga (y, g) being independent of g, a property that in two dimensions 
is obtained for v = Z. In such a case one obtains a scattering angle given by 



(39) 


which leads to the following expression for cr{g,x) 





(40) 


With this result, the integral / yields 


I = - 


/2k TT^n^ 

V m 2/3^ 


(41) 


where F is given in by Eq. (j64p in the Appendix. Thus the heat flux and thermal 
conductivity are given by 



(42) 


and 



(43) 


respectively, with F r; 0.212383. It should de emphasized that for this model, the 
thermal conductivity grows linearly with the temperature, which is analogous to 
result in the 3D case. 

5 Discussion and final remarks 

In this work, the standard procedure to establish transport coefficients has been 
used in order to obtain the thermal conductivity for a dilute bidimensional gas. As 
mentioned above, the temperature dependence and the proportionality con¬ 
stant for the hard disks case are consistent with the result in Ref. [T], which is 
widely used as the low density limit when considering dense gases. 

The main contribution of this work, the thermal conductivity given by Eqs. ra 
and (I32]), allows the introduction of any molecular model in order to obtain numer¬ 
ical approximations for this coefficient in bidimensional dilute gases. Additionally, 
the thermal conductivity is here calculated for Maxwellian molecules where the tem¬ 
perature dependence coincides with the three dimensional case. It is of the authors’ 
opinion that the explicit and systematic fashion in which this work is presented 
will be useful for generalizations to other dimensions, intermolecular potentials and 
transport coefhcients. To the authors’ knowledge, this type of calculation has not 
been published in the available literature and could also be used as a baseline for the 
calculation of bidimensional transport coefficients for relativistic dilute gases which 
will be reported elsewhere. 
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Appendix 


In this Appendix, the integral / given in Eq. (1^ , is calculated following the proce¬ 
dure described in Ref. m- Firstly, a change of variables is introduced (to relative 
and center of mass velocities): 


g = c — Cl G = Cl + c, 


such that 


G + g 


Cl = 


G-g 
2 ’ 


(44) 


(45) 


and the corresponding Jacobian is J = 2 By introducing these new variables in 
Eq. (1^ . I can be expressed as 


I = 


35 


j-(o) jio) 


( 


G^ + 2G-g + g 


0 


2kT 


- 2 


(c + g^ G-{]f^-gg) gad^ifGdPg. 


(46) 

Introducing the parameter /3 = m/2kT , and the equilibrium distribution function 

(V) 


one has 


M(?/3)"/exp 


G +2G-g+g^ 


liG^+9^)] K 

G-ig'g'-gg)] ga {x,g)dx<PGcPg. 


) “ 2 -I- 


(48) 


The integration with respect to G can be separated into two terms as follows 
{G^ + 2G-g + g^ 


exp 


~2^ 


= / exp 


-I- / exp 


2^ 


2^ 


P 


P 


- 2 


G^ + 2G-g + g^ 


G^ + 2G-g + g^ 


(g + pY G- - gg) ga (y, g) dxd^Gd^g 

G ■ (sY - M) 50 - ix, g) dxd^Gd^g (49) 
G ■ (^(f - 55)1 5CT (x, g) dxd^Gd^g. 


- 2 


G- 


- 2 


9- 


For the first term (second line in Eq. (HlUl '). since 


G- 


G-{^!f -gg) =Y1 GiGj {g[g'j - gtgj) 


8 





































one has 


idWj - 9i9j) / exp 


* j=i 


"2^ 


G2 + 2G-g + g2 


- 2 


G^GjtfG 


= Y1 (sWj - 9^9j) J exp 

jj'=i 


P f-i2 
~2^ 


(G^+9^) 

P- -7^ - 2 


which vanishes since the integral is proportional to 5ij and 

{9W] - 9 ^ 9 j) Szj = 5 '^ - 5 ^ = 0 . 

For the second term (third line in Eq. (Hm) !. considering that 


GiGjifG, (50) 


(51) 


9- 


G ■ -gg) =Y 1 9 iGj - 3*5^) 

jj'=i 


the integration with respect to G yields 
2 „ 


Y1 9^ iaWj - 9^9j) 9 J exp 

jj=i 

= ^ 9i {gWj - 9^9i) 9 / exp 




i,l=l 


= ^^9^ {gWj - 9^9J) 


P (-i2 
"2^ 


P 


G^ + 2G-g + g^ 


G^ + 2G-g + g^ 


- 2 


- 2 


Gjd^G 


Gjd‘\52) 


P 


and thus 


9o 


i,j=i 


I = 


rP/3 


^ (“ 2 ^^) 9 ^ 939 ^ ix,g) dxd^g, (53) 

i ,7 = l ^ 2 


The integration over the scattering angle y can be performed by aligning the y axis 
with the initial relative velocity such that 


and 

Since 

we have 


5 = 5 ( 0 ,1), (54) 

^ = 5 (sinx,cosx). (55) 

p‘2'k p'2'k 

/ o- (5, x) [gWj - 5 » 5 j] gidjdx = -9^ ^ (5, x) sin^ xdx ( 56 ) 

Jo Jo 




which is Eq. m- 


9 
























A specific molecular interaction model needs to be introduced in order to perform 
the X integration. For hard disks o’ = f |sin and thus 


I = - 


d f ( P 2\ 5 


-(f) 


sin^ Xdxd 9 


(58) 


Integrating over the scattering angle |sin (x/2)| sin^ x^X = 32/15, one obtains 




(59) 


which is the expression used in the main text to obtain the thermal conductivity for 
this model. 

On the other hand, for Maxwellian molecules, since x is given by 


2 ds 


X = 7r- 


(60) 


(l - s2 - (^) 

mg^ \s/ 

where s = b/r and Smax is the positive root of 1 — — 2$ {b/s) Img^ = 0, one has 


2 ds 


X = TT ■ 






= TT I 1- 

a 


(61) 


where = 1 + (2K/m) (l/g^b^). Equation (1611) is the expression for x given in Eq. 
dMl)- Solving for the impact parameter one obtains 


b{x) = 


2k 

mg'^ 


- 1 


- 1/2 


(62) 


.7^ - X. 

which leads to Eq. (HHll after differentiation. Introducing a {g, x) in Eq. (nTTll yields 

n^TT^ /2 k, 


I = - 


2/3^ V TO ’ 


(63) 


where we have defined 


r = y sin2x(x(27r-x)) dx- (64) 

The integral in equation (|M1) can be solved numerically to obtain F Ri 0.212383. 
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